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Some explicit Krein representations of certain 
subordinators, including the Gamma process 

C. Donati-Martin and M. Yor * 


Abstract 

We give a representation of the Gamma subordinator as a Krein 
functional of Brownian motion, using the known representations for 
stable subordinators and Esscher transforms. In particular, we have 
obtained Krein representations of the subordinators which govern the 
two parameter Poisson-Dirichlet family of distributions m- 

Mathematics Subject Classification (2000): 60G51, 60J55 

1 Introduction 

(l.a) The aim of this paper is to represent explicitly a particular class of sub¬ 
ordinators (S'tU > 0), i.e. increasing M+-valued Levy processes, as inverse 
local times of M+-valued diffusions. This problem was raised in Ito-Mc Kean 
|Sj, and completely solved, in a theoretical manner, thanks to Krein’s repre¬ 
sentation of strings, by Knight cn and Kotani-Watanabe m independently 
and simultaneously in 1981-1982. 

Roughly, if the Levy measure z/ of {St] t > 0) admits a density h with respect 
to the Lebesgue measure: ^{dy) = h{y)dy and if h{y) = y{dx) exp{—yx) 
for some positive measure p, then the above mentioned problem, which we 
shall call the Krein representation problem, is solved in the affirmative. See, 
besides uni and nn, a number of other papers dealing with this question, 
e.g. Bertoin (in. 0). Kiichler ^21, Kiichler and Salminen PI- 
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(l.b) In this paper, rather than discussing this problem in general, we solve 
it for the two parameter family of subordinators (5'“’^; t > 0) whose Levy 
measures are 

= C'—0<a<l,/?>0 
(for a = 0, /3 is assumed > 0). 

The constant C plays a simple role and might, a priori, be suppressed from 
our discussion; however, it is in fact very helpful to keep this further degree 
of freedom in order to consider the most convenient local time at 0 (for 
the underlying diffusion), which, as is well known, may be chosen up to a 
multiplicative constant. 

(l.c) It goes back at least to Molchanov-Ostrovski I2DI that for 0 < a < 1, 
the stable subordinator of index a, may be realized as the inverse local 
time of a Bessel process (which we shall denote as BES{—a)) of dimen¬ 
sion 6 = 2(1 — a). Since, by inspection of the Levy measures: Va^p^dy) = 
exp(—/3 ?/)z/q^o(<^I/), may be obtained as an Esscher transform of it 
is natural to look for a Girsanov transform of BES{—a) whose inverse local 
time is distributed as 5'“’^. As we show below, this is indeed the case with 
BES{—a,P J,) the downwards BES{—a) process with ’’drift” /3, following 
the terminology of Watanabe HQ, inzi and Pitman-Yor IZQ. 

Finally, the case a = 0, /? > 0, which corresponds to the Gamma process 
(*5'°’^; t > 0) is obtained from the case a > 0 by letting a —>• 0 in a suitable 
manner. 

(l.d) It is also natural to look for some representation of these subordina¬ 
tors as continuous additive functionals of a Brownian motion (i?^; s > 0) 
taken at the inverse local time (r/"; t > 0) of that Brownian motion. This is 
done expressing the Bessel processes with drift in terms of their Feller rep¬ 
resentations, using scale functions and speed measures. For example, it was 
remarked in Biane-Yor j3] that: 

(1.1) ( r du, t > 0) “=’ > 0), 

Jo 

and that all symmetric stable Levy processes may be obtained in this man¬ 
ner, replacing the even power by ’’symmetric” powers i.e.: a{x) = 

sgn{x)\x\^~‘^. For example, (y/J"* t > 0) is a standard Gauchy process. 
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More generally, every asymmetric stable Levy process may be represented in 
a similar manner from Brownian motion. Altliongli in the seqnel, we shall 
also extend (HU to present in general (5'“’^; t > 0) in terms of Brownian ad¬ 
ditive functionals taken at (r*;t > 0), these presentations are not so simple, 
and we prefer to those the Krein representations evoked in (l.c). 


(l.e) The rest of the paper is organized as follows. Our main Krein rep¬ 
resentation results are presented in Section 2. The full proofs are given in 
Section 3. The Brownian additive functional representations are discussed 
in Section 4. Finally, in Section 5, we also give some Krein representations 
of the symmetric Levy processes on R (without Gaussian component) whose 
Levy measures are given by: 




exp{-p\y\) 


(0 < a < 2, /3 > 0). 


(l.f) As an end to this introduction, let us point out that, the subordinators 
we represent here, and their symmetric counterparts, are arguably the 
most studied and used among Levy processes, and this, for the following rea¬ 
sons: for a > 0, is obtained by Esscher transform (see [71I2SI) from the 
fundamental stable (a) subordinator, hence it ’’retains” some scaling prop¬ 
erty, while the Gamma process (i, HZi, EH, EH, m) and the variance- 
gamma processes (HZl. Cn], HH) have some fundamental quasi-invariance 
properties, which make them comparable, in some respect, to Brownian mo¬ 
tion with drift. 


(l.g) Finally, we refer the reader to a systematic compendium |H of constants 
Ca related to various choices found in the literature of local times with respect 
to BES{—a). We also intend in |1] to study a number of properties of the 
symmetrized BES{a, [3 1) processes with a view towards a discussion of Krein 
representations for the variance-gamma processes. 


2 Definition of the BES{—a, 1) processes and 
main Krein representation results 

(2.1) We hrst recall that for 0 < a < 1, BES{—a), which we call the Bessel 
process of index —a, (or dimension d = 2(1 — a)) is the R+-valued diffusion 
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with infinitesimal generator 


1 (P 1 — 2a d 

2 dx‘^ 2x dx^ 

which is instantaneously refiecting at 0. We denote by Pi its distribution 
on C(M+,M+), where Rt{oj) = a;(t), IZt = cr{Rs] s <t} and Rq = x. 

For the sequel, it will be convenient to introduce the parameter 6 = for 
/3 > 0. We may now define BES{-a, P1) as the diffusion with law (Pi“"^’^'^), 
obtained by Girsanov^ transform from BES{—a): 

(2.1) = l^^exp(C„(2/3)% - 

Ra{0x) 

where Ka{,x) = x^Ka{x),x > 0, which satisfies Pa(0) = 2 ““W(q;) and 
{h',t > 0) is a choice of the local time at level 0 of P made so that the 
following holds: 

Ito’s excursion measure n^{de) associated with BES{—a, (31) together with 
this choice of local time may be described as follows: 

a) n^{y{e) G dv) = ^ exp(—where V{e) denotes the lifetime 

of the generic excursion e. 

b) Conditionally on 1/ = n, the process {e{u),u < v) is distributed as a 
Bessel bridge of index a, and length v. 

Note that this description is valid, in particular, for (3 = 0. 

The conditioned diffusions BES{—a, (31), started at a: > 0 and killed at Tq 
are described in [21]. See also Watanabe |SI1E2103|. 

We now discuss some immediate consequences of these choices: from (a), 
it follows that: 

Pj““^’^^(exp(-Ar„)) = exp (^u2^V{a + 1) ^ ^(1 “ 

and elementary computations yield: 

P^""^’^^(exp(-Ar„)) = exp {^u2°‘T{a + ^^ {(A + /?)“ - 

= exp(-u2‘^ . ^ J (A + /3)° 

\ sm(7rQ;) / 

^We shall also use the well-known general fact that such an absolute continuity rela¬ 
tionship extends with t replaced by any stopping time T on the set (T < oo). 
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As a consequence of these computations, it follows that, on one hand: 


( 2 . 2 ) 


"^(exp(-Ar„)) = exp (-u2°' . X 

\ sin(7rQ;) 


and from (EH) taken at x = 0 and t = Tu, we can determine the constant Ca'- 


a 


71 

sin(7ra;) 


On the other hand, 

lim h;J""^’^^(exp(-Ar^.)) 

a. -> 0 

(2.3) 


exp (-n{ln(A + /?) - ln(/d)}) 
1 

(1 + A//?)“ 


Thus, assuming that we may represent the LHS of as EQ’^^(exp(-Ar„)) 
for the law Pq’^^ of some diffusion BES{0, f3[), instantaneously reflecting at 
0, we will thus have obtained a Krein representation of the Gamma process 
with parameter /3, i.e. the subordinator whose Laplace transform in A, at 
time u, is given by (Q. 

Theorem 2.1 1) For every (3 > 0, there exists a diffusion BES{0, P 1) on 
M+, which is instantaneously reflecting at 0, and whose infinitesimal gener¬ 
ator on (0, cx)) is given by: 


1 d^ 

2 dx^ 





d 

dx 


2) For some choice of the local time at 0 of BES{0, fl 1), the corresponding 
ltd measure of excursions may he described as: 

a) n^{y{e) G dv) = ^e'xp{—j3v)dv 

h) Conditionally on V = v, the process {e{u),u < v) is distributed as a 
Bessel bridge of dimension 2, and length v. 


3) Let X > 0, the following relation holds: 


(2.4) 
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Note: The following relation, where we have extended {1/Kq) to [0, cxo[, with 
(l/it"o)(0) = 0, is equivalent to (12.411 : 


(2.5) exp(/3()Pi">’‘’‘k.. 

-fvo(.[7tttAToj 

The following theorem exhibits an absolute continuity relationship between 
the laws and from which several important Laplace trans¬ 

forms may be immediately obtained. 


Theorem 2.2 Let /3,7 > 0, and x > 0; then, there is the relationship: 

( 2 . 6 ) 

p(o),( 0 + 7 )l|„ ^ Ko(,/2{P +'i)R,) K„{,/2px) —/ -.pdIciI 

K„{VWR,) K„(Pp3TT)x) L ^ ' “ 


/3i| 


|7^t• 


In particular, for a; = 0 and t being replaced by ti, formula (ESI) simplifies 


as: 

(2.7) P„ 


and 

( 2 . 8 ) 


(0) 


.(/3+7)t| i^o(\/2(/3 + 7)i?t) 


= (1 + ^)kxp(-7rOPf 


Note, also from (EH), that the measure: 


K„{x/pR,f 

does not depend on jS. 

Here are some consequences of (EH, (EH, (EH: 

- As a consequence of (EH, one obtains the Laplace transform of 
To = inf{t; Rt = 0} under 


(2.9) 


pW’^kexp(-7To)) 


Ko{^/2{(I + 'y)x) 

KoiVWx) 
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In fact, Tq is distributed^ as a GIG{0;x, s/^) variable, i.e. with den¬ 
sity: 


1 

2KoixV^) 


exp( 


1 

2 


-|- 2l3t))dt. 


As a consequence of dZHD, we obtain our main result: 

Corollary 2.3 (rp / > 0) is, under , a gamma process with parameter 
[3, i.e. : I >0) is a standard Gamma process: 


( 2 - 10 ) Ej°^’^^(exp(- 7 r 0 ) = ^ ^ for 7 > 0 . 


3 Esscher and Girsanov transforms 


This section is devoted to the proofs of the assertions contained in Section 2. 

(3.1) To begin with, we explain how, once we have obtained a Krein repre¬ 
sentation for a subordinator {St), we can obtain a related one for a second 
subordinator {St) whose law is an Esscher transform of that of {St). We 
consider a positive diffusion {Xp, t > 0) with 0 as an instantaneous reflecting 
boundary. We denote by Lt a choice of its local time at 0 and {ti : I > 0) 
the corresponding inverse local time. We are looking for a diffusion X on M_|_ 
with laws Px such that: 


(3.1) 


Po\tt, = exp(t/i(o)t( - ar„)Folx^„ 


for a > 0, where tjj denotes the Laplace exponent of the subordinator (r„; m > 
0). From ()3.1|1 . the inverse local time r of X is the Esscher transform of 
T dehned by: 


(3.2) 


Eo[exp(-Ar^“))] 


Eo[exp(-(A a)Tu)] 
Eo[exp{-aTu)] 


The Levy measure m^°‘\dx) of is related to the Levy measure m{dx) of 
T by m^°'\dx) = exp{—ax)m{dx) (see j7], [211, Chapter VII, Section 3c). 

To dehne Px, we need to compute the martingale 


Ms = E[exp{-aTt)\Ps], s < n 

^For many references to, and applications of, GIG variables, see Matsumoto-Yor m 
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Ms = exp(-as)Ex4<^) [exp(-art_L4^))] 

Now, we have: 

[exp(-ar„)] = [exp(-aTo)] Eq [exp(-ar„)] 

and therefore: 

TT = ^oi{Xs) exp(^/>(a)L, - as) 

Mq 

where 994 denotes the function dehned by 

(3.3) ^aiix) = E^[exp{-aTo{X))] (<^ 4 ( 0 ) = 1) 

where To{X) is the hrst hitting time of 0 by X. More generally, for a; > 0, 
we dehne the law Px via the absolute continuity relation: 

(3.4) Px\rt = exp(^(a)Lf - at)P^|^,. 

The generators L and L of X and X respectively are linked by: 

1 = L + ^^—. 

Lfai (x) dx' 


(3.2) Examples: 


1. Px = Px 0 < q; < 1. Then, 

(Pai{x) = CaKaiV^x), ^/’(a) = . J {2a)°‘ 

sin(7rQ;) 

and Px = 

2 . Px = Pi““^’^\ 0 < a < 1, /3 > 0. Then, 

vw = -^{(2o+ 2pr - (2/3)“} 

K^{^x) sin(7ra) 

and Px = pi-“)’(d+a)i_ 
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From [H] and I2DI. we know that the inverse local time of BES{—a) is a stable 
snbordinator of index a. It follows from snbsection (3.1) and Example 1 
above that the inverse local time of BES{—a,P |) is a snbordinator 5"’^. 
The description of Ito’s excnrsion measnre of BES{—a,f3 J.) follows from 
the description of the excnrsion measnre for BES{—a) given in and the 
Esscher transform. Note that we have not chosen the same normalisation for 
the local time as in [3]; we have the following relation between the two local 
time processes (from ()2.2j) and snbsection 3.2 in j^): 

by _ ( 2 ° T ( a ))^ by 

* “ 2(1-a) ‘ 

where denotes the local time considered in IHj. Note that there is a 
mistake in jS] after (3.i) dne to the identihcation of (the inverse local 

time of the reflected Brownian motion) with r* (for the Brownian motion) 
instead of Tt /2 and the correct formnla (see p.45, after (3.i)) is 


(3.3) Proof of Theorems 12.11 and 12.21 

1) We consider the diffnsion X on [0, cxo[ with generator 


1 


K 


nXX + FZ + -JZ 


d 


2 dx^ ' \2x ' '' ' iFo ^^ J dx 
A pair {s{x),m{dx)) of scale fnnction and speed measnre is given by: 

m{dx) = 2xKl{\/^x)dx. 

1 


s X = 


/o yKl{^y)' 

Note that KQ^y) ~o ln(2/|/), thns 


is integrable in 0 and s(0) = 0. 


yKiiV^y) 

In order to obtain the behavior of X at the bonndary point 0, we shall apply 
the criterinm of Rogers-Williams 123 V.51] to the diffnsion in natnral scale 
Y = s(X). The speed measnre of Y is given by: 


niY^dx) = 


(s' o s“^(x))- 


-dx. 


Then, xniY^dx) < oo implying that To{Y) < oo Px p.s. for all x > 0, 
and fg_^m(dx) < oo ensnring that 0 is reflecting. 
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The description of the ltd measure follows from that of BES{—a, /3 J,) letting 

a —>• 0 . 

The absolute continuity relation Q is given in Pitman-Yor [25. We can 
also obtain it from m and: 


Thus, 


= (^) ^> 0 - 

\ X j K^{ex) 


and letting a —0, 

The absolute continuity relation (EH) follows from the example (3.2) 2. with 
a = 0. In this case, 

/ N „-Yo(\/2(/d -1- a)x) _ a. 

iPai{x) = C ^a) = lnl + -. 

Kq[^/2I3x) P 

The other formulas follow easily. □ 


4 Brownian representations of the subordi- 
nators 5“’^ 

We keep the same normalisation as in Section 2, i.e. the Laplace transform 
of a stable process is given by (Q 

Proposition 4.1 1. The stable subordinator S°‘ can be represented by the 

Brownian additive functional (see m 


Bu \“ ^ du, s > 0, 
and Tt is the inverse local time of B. 


where 

(4.1) kl„(5) := r 

Jo 
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2. The subordinator 5'“’^ can he represented as: 


(4.2) 


hc,{\Br\) dr. 


where ha{x) = 2{s^^{2x)fK^{s^^{2x)) and 


(4.3) 


SniX = 


dy 


lo yKl{y) 

3. The Gamma process 5'°’^ has the following representation: 

(4.4) h{\pu\) du; t >0 
with 

(4.5) h{x) = 2{G-\2x)fK^{G-\2x)), 

where G~^ is the inverse of the increasing function G given by: 

(4.6) G(z) = 


dx 


xKq{x) 

Sketch of Proof: 

l)The first point has been obtained in [21. More precisely, 


(4.7) 

where 


E[exp{-^Aa{Tt))] = exp(-tCaA“), A > 0, 


TT 


Crv 


a 


asin(7ra) \r(Q;) 

2) By a Girsanov transform, we can find a diffusion X such that the subor¬ 


dinator 


: = 


^Tt(X) 


IX 


I--2 


du 


is a Esscher transform of the stable subordinator Aa{Tt) and thus is dis¬ 
tributed as 5'“’^(up to a constant). Then, we write the diffusion in natural 
scale Sq,(X) as a time change Brownian motion. We can also start from 
BES{a,l I) and use a time change method. 
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3) By an application of Ito’s formnla, we can prove that: 

G-\m\) = Y(^j\h{m)ds^ 

where F is a BES{0, 1|) process from which we can dednce that Tt{Y) 

2h(\Bs\)ds where Tt{Y) is the inverse diffnsion local time of F, dehned 
with the speed measnre rriY^dy) = 2yKl{ii)dy. 


5 Some complements 


5.1 The case o: = 1/2 

Let ns go back to the representation (Q for the snbordinator with Levy 
measnre fc^exp(—In the case a = 1/2, (14.21) becomes: 


(5.1) 


dr 


Iq (l + 2|/?r|)2 


We now explain how this representation may be rednced to the more classical 
one, as given in Q below. For this pnrpose, we recall that the stable 
snbordinator of index 1/2 can also be realized as the distribntion of {Ty, f > 0) 
where Tt denotes the first hitting time of t by a Brownian motion, starting 
from 0. We have the same interpretation for its Esscher transform. Indeed, 


ln(l + 2|i?t|) 


G sgn{Bs)dBs ds ' 

[ {i+m\) Jo ( 1 + 215 , 1 ) 2 + 



ds 


(1 + 2|5, 


+ Lt 


where > 0) is a Brownian motion with drift 1. 

From Skorokhod’s lemma, we derive: 


Lt = snp < 


dr 


lo (1 + 2\ Br\y 


Therefore, we obtain: 

r' dr 


(5.2) 


'0 (l + 2|Br|)‘ 
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Thus, from (I5.1|l . / > 0) is a subordinator with Levy measure C exp(— 

and Laplace transform 


y\ d.y 



Of course, the above Laplace transform could also have been computed from 
the Laplace transform of Ti and the Cameron-Martin formula. 

5.2 Symmetric Levy processes 

From the Brownian representation of the Gamma process, we can give a 
Brownian representation for the symmetric Gamma process (or variance 
gamma process) distributed as (71 (t) — 72 (^)G > 0 ) where 71 and 72 are 
two independent gamma processes. More generally, we have: 

Proposition 5.1 Let Si and S 2 be two independent subordinators with Brow¬ 
nian representations: 



Then, we have the Brownian representation for (S'i(f) — S 2 (t);t > 0) as 
follows: 


(5.3) 




Proof: It suffices to write 



dsip{Bs) = / ds\Bs>o)Ti{Bs) - / dsl{^Bs<o)T2{-Bs). 



Then, we use the representation 


(5.4) 
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where and are (as a consequence of Knight’s theorem) two indepen¬ 
dent Brownian motions^ ; see, e.g., proofs of the arc sine law for Brownian 
motion, inspired from D. Williams jSl] (see i, mi). 

Then, we can write 


r dsip(B,)= T" 

Jo Jo Jo 

where 

rrt pTt 

= / dsl(^Bs>0)ds-, = / dsl(Bs<0)ds, 

Jo Jo 

and moreover, from D, we also learn that: 

-4^.+^ = 'rt/2(/3(+)), = Tt/2{f3^~^). 

Finally, we have obtained the result. □ 


5.3 The Ito measure of BES{—a,Pl) 

The description of the ltd measure given in Section 2 relies upon one of the 
descriptions of the ltd measure ria of BES{—a) given in |3]. There is a 
second description of n conditionally to the maximum of the excursion. 

Proposition 5.2 a) Under n^, the distribution of M satisfies: 


(5.5) 


nd^{M >x)= 2(2/3)“ 


Ka{y/^x) 

laiVWx) 


b) Conditionally on M = x, the maximum is attained at a unique time R and 
the processes {et;0 <t < R) and {ey-t^ 0 <t <V — R) are two independent 
BES{a, processes, starting from 0, stopped at their first hitting time E 
ofx. 

Proof: a) From the description of given in section 2, we have: 


n 




Ifitr-)] = 2“r(a +1) /” lf{Ar)]dv 


^However, we emphasize that, knowing B, only the reflected Brownian motions 
and I are accessible, and not and /jK). 
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where is the distribution of a Bessel bridge of index a and length v. 

"?[/(A/")l = 2"r(a + 1) r 

Jo ^ 

Now, from j22l Theorem 3.1] 

ni“)[0(r)]=c„E[0(.R)(M)-2“] 


where the process R is dehned by 


Rt — {T + T) 0 < t < 1 


and Y is the process connecting the paths of two independent BES{2 + 2a) 
processes R on [0,T] (hrst hitting time of 1) and R on [0,T] back to back, 

i.e. 


Yt 


Rt t <T 
Rt+t-i ^ < t < T + T 


Then, 


M := sup.R, = (T + f )-^/2 

S<1 


and Ca = 2“r(Q! + 1). 
It follows that: 


nS“)[/(nM^)]=c«E 


/ 


T + T 


{T + ry 


Thus, 


na[/(M^)] = 2“r(« + l)c„E 


= 2 "r(a + l)co 


exp(-/3M;(T + f)) 


w 


a+l 


f{w)dw 


f{w) 


w 

coo 


a+l 


E 


exp{-(3w{T + f)) 


dw 


[2"r(a + l)]-^c^ r ^ ] dw 


( 2 /d)“ 

2{2(3y 


fiw) 


'o wlaiVWw) 

r f{y^) 

Jo yiliVWy) 


dw 

dy 
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Therefore, we obtain 


hence: 


n^a{M > x) 


2(2/3)" 

2(2/3)" 


1 


Jx yiaiVWy) 

K^iVWx) 

laiVWx) 


dy 


We also refer to 123 for related computations. 
As a verihcation, we can let (3 —>■ 0 to obtain: 


na{M >x) = 22"r(«)r(« + 1)^. 

This agrees with the description of the distribution of M under given in 
Biane-Yor jS]: 

h„(M > ^) = ^ 

and from our normalisation given in Section 2, 


Ua = 2^"r(a)r(a + l)ha. 


b) We refer to [21] for the dehnition of Bessel processes with drift BES{a, /3|). 
The description of conditionally to M, follows from the description of 
Ua (see 121), the relation n'^{de) = exp{—pV{e))na{de) and 

= C'«x-"4(0a;)exp(-/3Tjp("V.. 

(see Proposition 3.1 in mi). □ 
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